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Abstract 

In a recent study by Dittmar and Dreiner it was shown that with appropriate 
selection cuts the signature of events containing two charged leptons and 
missing energy represents the best chance of detecting the Standard Model 
Higgs scalar in the mass range between 155 and 180 GeV, the primary 
decay of the Higgs being into pairs of charged gauge bosons. The largest 
background to this channel is due to irreducible W^W~X production. In 
the present paper we calculate the contribution of events of the type bg — >• 
tW^ — i> bW^W~ — > b[^[ ^vivni, which have not been considered yet within 
the new selection strategy. We show that the yield of this background is 
rather large, at the level of that produced by W~^W~, tt or tbW^ events and 
thus needs to be incorporated in future experimental analyses. However, 
we find that its inclusion will not spoil the possibihties of Higgs detection 
in the above mentioned channel at the Large Hadron Collider. 
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1. Introduction and motivations 



In a recent paper by Dittmar and Dreiner lU (see also Ref. 0) it was pointed out that 
the signature of events with two charged leptons and missing energy /momentum at 
the Large Hadron Colhder (LHC) represents the best chance of detecting the Standard 
Model (SM) Higgs boson in the mass range 155 GeV 1^ Mh ^ 180 GeV. In Ref. 0, 
simple selection criteria were outlined, which should allow one to extract the Higgs 
decay channel H —>■ W~^W~ from the non- resonant W~^W~X production (where X 
represents possible additional particles in the final state) with a signal-to-background 
ratio of about one-to-one, thus allowing a 5^10o" detection with only 5 inverse picobarns 
of integrated luminosity C = J Ldt. The appealing prospect is that this significance 
can be achieved in less than one year of running of the CERN machine at the initial 
low luminosity L = 10^^ cm~^ sec~^. Indeed, this is a clear improvement compared 
to the Higgs search strategy based on the decay mode H — > ZZ* — >• four charged 
leptons, which was the detection channel exploited even in the most recent experimental 
simulations 0, ^ for the mentioned Higgs mass range. This is evident if one considers 
that in order to disentangle a 5a signal in the latter case at least 100 fb~^ are required. 

The first studies of the H W^W~ decay mode Q in the context of Higgs 
searches at the LHC date back to Ref. and to the 1990 Workshop [§] for a LHC 
with a/s = 16 TeV. Further analyses were subsequently performed, in Ref. 0. In var- 
ious instances, also several signal-to-background studies were carried out (see Section 
2 of Ref. |I| for a review). The unanimous conclusion was that the H —>■ W^W~ —>■ 
i~^i~h'ih'i/ channel (with = e,/i) should provide a useful tool to detect the Higgs 
boson in the mentioned mass range, though more appropriate analyses (including hadro- 
nisation and detector effects) were recognized to be needed to support those (mostly 
parton level) results. 

This was done in Ref. |^, by using the Monte Carlo (MC) program PYTHIA ||10|| . 



Further refinements were also introduced there, which were not included in the previous 
literature. Namely, (i) the inclusion of W"^ — > r^z/^ i'^uiUj. decays (with i = 
e,/i); (ii) the simulation of the background due to gg tbW^ events [|ll], |12|; (iii) 
cuts previously employed 0, |^ were further supported by new constraints, introduced 
mainly in order to discriminate against the 'irreducible' background from continuum 
production of W~^W~X events. 

It is the purpose of this letter to provide additional material to motivate the ex- 
ploitation of the H W~^W~ channel in Higgs searches at the LHC, as we have studied 
the irreducible background due to 

bg ^ tW^ bW+W bi+l'-ueUi, © C.C., (1) 

'single-top' events via fef^-fusion (also called 'tVT^-production'), which was not consid- 
ered in Ref. |l|, and we will show that this can be reduced to a manageable level by 
the same selection criteria recommended in |l|]. In fact, for completeness, we have also 
computed the yield of the process 

bg bW+W- bi+i'-UfUfj ® C.C., (2) 

involving all the tree-level graphs producing the final state bi'^ i' ^ uiup : that is, not only 



the single-top ones isolated in reaction (|^) but also all the other diagrams contributing 
at the perturbative order 0{aj^as)^. 

With respect to the analysis performed in Ref . , we will adopt two simplifications, 
which we believe will not spoil the validity of our conclusions. First, although we will 
implement the same cuts considered in Ref. [|I|, we will confine ourself to the parton 
level only. However, since at lowest order the final states of reactions ([l|)-(|^) involve 
only one hadronic system (i.e., the 6-quark fragmenting into hadrons) whereas the 
Higgs signal H — > W~^W~ — > ~z/^z/£/ is purely leptonic, we expect the effects of 
hadronisation not to modify drastically the parton level dynamics. Second, we will 
only discuss the channels W^W~ ~z/£Z/^/ with £, = e,/i, thus neglecting the 

case of ly^-decays into tau leptons via the three-body channels r^Vr C^viVr. 

This is done to simplify the description at parton level (especially in the case of the 
complete process (|])), as in this way we can avoid to calculate complicated two-to- 
seven and two-to-nine body subprocesses. In practice, contributions involving r-decays 
amount to ~ 1.9% of the total ~ 7% leptonic branching ratio of """ly "-pairs, so that 
the bulk of the produced W^W~ events are indeed included in our study. In general, we 
stress that we are here only interested in the relative rates of signal and background and 
we expect that the implementation of a full Monte Carlo simulation and the inclusion 
of the — >■ TUr decays will presumably affect both in a rather similar manner. 

The reason for studying processes (|l|)-(@) as a potential background in Higgs 
searches in the two leptons plus missing energy channel is that single-top production via 
process (|l]) has very large event rates at the LHC, as its total cross section amounts to 
55 — 60 pb at -y/s = 14 TeV (see later on), thus being cornparable to that of the process 
gg — ^ tbW^ considered in Ref. [||] (see, e.g., Ref. [|1^, |13[p- Furthermore, we stress that 
compared to the final state tbW^, which eventually yields the signature bbW~^W~ , that 
of reaction (|^) (and, more generally, of the complete process (^) can boast only one 
additional particle with respect to the Higgs signature (this rendering its reduction less 
effective than that of tbW^ events, which have two additional jetsi). In fact, the latter 
is produced at lowest order via gluon-gluon fusion into an on-shell Higgs boson, through 
a top quark loop |M: gg ^ H ^ W~^W~ —>■ £^£ 'fiVi'. However, we notice that the 



ii'-factor of Higgs production via gg-insiori has been shown to be very large, around 

two |1|, ITeHg in the mass range 155 GeV ;£ ^ 180 GeV (and outside, as well 

In particular, a large component of the next-to-leading (NLO) order corrections to the 



gluon-gluon fusion mechanism of Higgs production is due to the real radiation |18[ of a 
quark or gluon, so that also signal events are naturally accompanied by an additional 
detectable jet inside the detectors. 

For reference, we recall that the matrix element of process (|l]) was already computed 



^The symbol © C.C. means that we have calculated also the charged conjugated processes initiated 
by 6g-scatterings and these are included in all event rates presented in the following Sections. 

•^Note that the leading order (LO) rates of the gg ^ H signal for 155 GeV ~ Mh ~ 180 GeV vary 
between 10 and 8 picobarns, approximately. 

*In this respect, we should mention that an extensive compilation and a detailed discussion of pro- 
cesses involving single-top production at hadron colliders has recently been given ||l2|. In particular, 
according to the classification of Ref. [Q, there are six of these different hard parton scatterings. How- 
ever, process is the only one contributing at lowest order to the irreducible background W^W~X 
with one additional particle in the final state (i.e., X = 6), as the others always produce a second 
(light) jet. 



in Ref. |T^ and first studied in the context of Higgs searclies (and of W^W~ physics, 
as well) in Ref. pO| (for its relevance in the case of top-quark physics, see Ref. [[12| ). 
However, only the invariant mass region Mtw = VS > 850 GeV was considered there, as 
a background to signatures of heavy Higgs bosons decaying into longitudinal polarised 
W+W--paiTs pTf. 



The plan of this paper is as follows. In the next Section we give some details of the 
calculation. Section 3 is devoted to a discussion of the results. Our conclusions are in 
Section 4. 



2. Calculation 

The tree-level Feynman 'topologies' that one needs for computing processes ([l|)-(@) are 
given in Fig. 1. Once all the internal propagator are correctly inserted one gets a total of 
43 Feynman graphs (the single-top diagrams pertaining to reaction (0) can be obtained 
from the topologies 2 and 3). To calculate the corresponding amplitude squared we 
have used MadGraph and HELAS [^. The integrations over the appropriate 



phase spaces have been performed by using VEGAS The codes produced have 

been carefully checked for gauge and BRS invariance. Furthermore, the total cross 
section for process (|l]), obtained by selecting the only two graphs with on-shell top 
production out of those displayed in Fig. 1, has been compared against the results 



given in Ref. [T^ for the Tevatron and in Ref. pO[ for the Superconducting Super 
Collider (SSC), with identical choice of parameters, cuts (where applied) and structure 
functions, and perfect agreement has been found. The signal rates have been computed 
by using the program already adopted in Refs. [E^, ETI]. However, contrary to the case 



of Ref. where NLO rates were used to calculate the Higgs production cross sections 
via gg-fusion, and in line with Refs. Jl], we have used here the LO results. This has 
been done for consistency, as one-loop calculations do not exist to date for processes 
(|l|) and (0). It is however important to point out that whereas the complete corrections 
to the main Higgs production process via gluon-gluon fusion are large and positive [0 



those to the single-top process (|l|) are expected to be much smaller |jT2 



The 6-quark in the initial state of reactions (|l|)-(0) has been treated as a constituent 
of the proton with the appropriate momentum fraction distribution fb/p{x, Q^), as given 
by the parton distribution functions (PDFs). As default set of the latter we have used 
MRS(A) ||2^. However, as the PDFs of the gluon inside the proton are not so well 
know at medium and small x and since those of 6-quarks suffer from potentially large 
(theoretical) uncertainties (see, e.g., Ref. [^), we have produced our results in the 
case of other 4 sets of recent NLO structure functions, which give excellent fits to a 
wide range of deep inelastic scattering data (including the latest measurements from 
the HERA ep collider) and to data on other hard scattering processes. These are the 
packages MRS(A', G, Rl, R2) RDI, M. The QCD strong coupling a, entering 



explicitly in the production cross sections and implicitly in the parton distributions 
has been evaluated using the CERNLIB package at the scale n = a/S (i.e., the CM 
energy at parton level). We will discuss the fi dependence of the LHC rates in the 
following Section. Indeed, one should remember that also the value of as associated 
with each parton set represents in principle a residual source of error in the predictions 



of the different PDFs. However, the value adopted in each set is chosen to match the 
data during the extraction, so that we do not expect to be a significant source of 
uncertainty. 

In the numerical calculations we have adopted the following values for the electro- 
magnetic coupling constant and the weak mixing angle: aem = 1/128 and siii^ 9w = 
0.2320. For the gauge boson masses and widths we have taken Mz = 91.1888 GeV, 
Tz = 2.5 GeV, Mw± = MzcosOw ~ 80 GeV and Tw± = 2.08 GeV, while for the 
top mass we have used = 175 GeV [Q. All other fermions have been considered 
massless, including the 6-quark. In particular, the choice = has been maintained 
also in the Yukawa couplings of the theory. In this way, no diagram involving radiation 
of Higgs bosons off the 6-lines has been included in process @. For simplicity, we have 
set the CKM matrix element of the top-bottom coupling equal to one. In this respect, 
we recall again Ref. [12|, where it was shown that off-diagonal CKM matrix element 
subprocesses are negligible at the Tevatron. We do expect the same to occur at LHC 
regimes. 

Finally, as total CM energy of the colliding beams at the LHC we have adopted the 
value ^/s = U TeV. 



3. Results 

Our results are presented in Tab. I and Figs. 2-6. Note that for the time being we 
assume that no 6-tagging identification is exploited in events of the type The 
integrated luminosity adopted throughout the paper will be 5 fb"\ 

3.1 Selection cuts 

As event selection procedure we have adopted the same one exploited in Ref. [jl[, to 
which we refer the reader for a detailed discussion concerning the meaning of the various 
cuts. We only tabulate these here, in order to introduce a notation that will be used 
in the remainder of this paper (note that the two leptons i and i' must be of opposite 
sign). Following the same numerical sequence as in we ask (at parton level): 

1. pff > 10 GeV, or pfp > 20 GeV, 6i i' > 10°, for the transverse momentum and 
the separation angle of the two leptons; 

2. |?7^'^'| < 2, for the pseudorapidity of the two leptons; 

3. Eb < 5 GeV if 9bi,be' < 20°, for the energy of the 6-quark and the separation angles 
between the leptons and the 6-quark; 

4. Mw < 80 GeV, for the dilepton mass; 

5. > 20 GeV, for the missing transverse momentum of the event; 



6. 



< 135°, for the angle between the two leptons in the plane transverse to the 
beam direction; 



7. ii \ri^\ < 2.4 then < 20 GeV, for the transverse momentum and the pseudora- 
pidity of the 6-quark; 

8. |cos^^| < 0.8, for the cosine of the dilepton system with respect to the beam 
direction; 

9. the enforcement 10° < < 45°, for the same angle defined in 6.; 

10. M > 140 GeV, for the estimated invariant mass of the W'^W' -system; 

11. < cos^ < 0.3, for the angle between the lepton with the largest transverse 
momentum, boosted to the dilepton rest frame, and the momentum vector of the 
dilepton system. 

A few comments are in order before proceeding further, concerning the application of 
cuts 3. and 7. As for cut 3., according to our parton level implementation, background 
events are rejected if the 6-quark is energetic and is found near one of the leptons. On 
the one hand, we certainly expect the 6-quark to be very fast. On the other hand, 
we do not see a priori any reasons why the quark and the leptons should be created 
in collinear configurations. This is in fact confirmed by the spectra given in Fig. 2a. 
However, things would look quite different at hadron level. In fact, the jet produced by 
the bottom quark would have a finite size and the hadrons produced in the showering 
would carry only a fraction of the original parton energy. Although we miss these two 
aspects, we stress that the two systematics errors we introduce with our treatment 
do work in opposite directions, so to counterbalance each other. Cut 7. will have 
no effect on our signal rates, as we are considering here neither r-decay modes nor 
initial state QCD radiation, whereas for the background it will act directly on the b- 
parton. This corresponds to an overestimate of the signal, while we believe that the 
accepted fraction of background events will be predicted accurately, as the efficiency 
in reconstructing the 6-momentum from the hadrons should be rather high because of 
the clean environment (the two leptons) in which the 6-quark fragments. Whichever is 
the interplay between parton and hadron level, is anyway clear that it is cut 7. that 
will introduce a strong reduction factor on the background, as the &-jet will be easily 
detectable in pseudorapidity and will also have a large transverse momentum! (see 
Fig. 2b). 

3.2 Theoretical error 

As first step of our analysis we have compared the production rates of process ([l|) 
and (0) and found that in Higgs searches (that is, for the selection cuts l.-ll.) the 
additional contributions from the non-top diagrams of Fig. 1 are negligible. Therefore, 
in the following we will neglect them. 

As one of the possible means of estimating the uncertainty of the theoretical predic- 
tions on the gluon distribution (and hence the 6-one) we have calculated the cross sec- 
tion for on-shell single-top production via the two-to-two body process bg —>■ tW^(BC.C. 



^Note that Fig. 2b has been plotted after having aheady implemented the constraints 1.-6., and so 
will be in all forthcoming Figures. 



for the mentioned five sets of PDFs. The spread around the value obtained from 
MRS(A) (the set that we will adopt as a default in the following) is between —9% 
(from MRS(G)) and +3% (from MRS(R2)). This will represent throughout the pa- 
per the conservative estimate at present time of the uncertainty on the 6(7-fusion cross 
section into single-top quarks due to the parton distributions. Note that the above 



values roughly compare to those identified (for the same sets) in Ref. ||2^ for the case 
of gg-iusion into an on-shell Higgs, so that this helps in this context in carrying out a 
consistent signal-to-background analysis. 

Finally, the factorisation scale dependence (which quantifies our ignorance of higher 
order corrections) of the background rates via process (|l|) has been estimated by varying 
H in the range < jj, < 2^/§ when calculating the total cross section. We notice 

that, using MRS(A), differences with respect to the rate at /x = are less than 0.1% 
at /i = \/S/2 and —3% at fx = 2\/S. We have verified that similar effects also occur 
when other PDFs are used. Such numbers are rather small and presumably comparable 
with the experimental uncertainties!. 

3.3 Kinematics and event rates 

One should expect the impact of the cuts 8. and 9. to be similar on both signal and 
background, as can be noticed from Figs. 3 and 4, respectively. In fact, the shapes of 
the corresponding distributions are almost identicalQ. 

Not even the invariant mass M of the reconstructed (from the lepton and the missing 
momenta) W^W~ -system is helpful to discriminate the signal from the background (see 
Fig. 5). In fact, the background spectrum is almost entirely beyond the minimum value 
of 140 GeV implied by cut number 10. The discrimination power of such constraint is 
thus very limited, if not self-defeating. 

The only cut among those introduced in Ref. IQ to reduce the irreducible W^W~X 
background from W~^W~, tt and tbW^ events which is also effective against bg —>■ tW^ 
is cut 11., as can be appreciated from Fig. 6. In fact, the two charged leptons from the 
background have a rather large angular spread, so that the maximum of the background 
distribution is located around the value 0.6. 

The accepted event rates, for both signal and background, for a selection of six Higgs 
masses, are presented in Tab. I {^/s = 14 TeV and £ = 5 fb~^). When comparing the 
numbers in Tab. I one should bear in mind that the background rates there should be 
added to those given in Tab. 2 of Ref. |]l|. This should however be done after treating 
all background sources in W^W~X events on the same footing (i.e., consistently at 
parton or, better, hadron level). This is beyond our intentions and capabilities, our aim 
here is to make the point that background events from process (|l]) are large compared 
to the signal, as they vary between 11% and 22% of the Higgs rates, depending on 
the mass of the scalar. Therefore, their effect in the signal-to-background significance 
is of the same order as that of any of the three processes pp W~^W~, pp ti 
and pp tbW^ studied in Ref. |l[], especially considering the fact that our parton 



^Note that the scale dependence of processes producmg the final state tW^X at the Tevatron has 
been studied in Ref. JT2[ , where variations between —14% and +20% were quoted, for fi spanning over 
the range between fi = mt/2 and 2mf. 

^Please notice the arrow in Fig. 4 to indicate the maximum value of the signal at cos0 w 1. 



level analysis overestimate the signal by a factor of two (compare the numbers in our 
Tab. I to those in Tabs. 1-2 of Ref. in response to the application of cut number 
7.), while more accurately predicting the background rates. Finally, one should notice 
the effectiveness of the selection strategy based on the cuts l.-ll. against events of the 
type (|1|), as the overall reduction factor on this background is above 1000 ! 

4. Summary and conclusions 

In this paper, we have studied the yield of process bg tW^ bW^W~ bi^ i' ~ ueui/ 
(where i, i' = e, /i) as 'irreducible' background to the gg ^ H ^ W^W" i'^i "i^ii^e' 
signature of the Standard Model Higgs boson, which has recently been claimed as the 
most viable channel to detect such a scalar in the mass range 155 GeV ~ ~ 180 
GeV at the Large Hadron Collider. Although we have confined ourselves to the parton 
level only, we believe we have performed a consistent signal-to-background analysis, 
exploiting the same event selection procedure advocated in literature. (In particular, 
the shape of the parton level distributions used to disentangle the signal from the 
irreducible W^W~X noise resembles very closely those previously obtained at hadron 
level). This has enabled us to assess that non-resonant W^W~X events due to single- 
top production via fogf-fusion are rather numerous, and comparable to the rates of 
any of the reactions gg,qq —>■ ti, gg — > tbW^ and gg,qq — >■ W~^W~ , which have in 
fact been shown to represent the largest components of the total background to the 
Higgs detection channel in two charged leptons and missing energy/momentum. In 
contrast, bg bW^W~ bi^i~h'eh'e' events not proceeding via single-top diagrams 
are negligible. Therefore, we think that the production process bg tW"^ that we 
have studied should be included in the experimental Monte Carlo simulations which 
will be used in order to confirm or disprove the existence of the Higgs scalar of the 
Standard Model in the above mass range at the CERN proton-proton collider. We 
believe this to be particularly important, as the discussed signature does not allow 
one to reconstruct the narrow Higgs resonance (because of the neutrinos escaping the 
detectors). In fact, the presence of the latter will be established by an 'event counting' 
operation over a rather broad region in mass, where the H — > W^W~ signal and the 
bg — > tW^ background have a very similar shape. For the purpose of aiding future 
analyses, we make available upon request the electronic version of the matrix element 
for bg bW+W~ bi+ i' - ueup . 

However, we would like to conclude this study by stressing that the inclusion of the 
single-top background in tW^ events will certainly not spoil the chances of detecting 
the Standard Model Higgs in the advocated decay channel, and that the exploitation of 
the two charged leptons and missing energy signature remains crucial in Higgs searches 
at hadron colliders. 
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Table Captions 

[I] The expected signal and background number of events for 5 fb~^ at the LHC after 
the apphcation of the selection criteria discussed in the text. Only the two-body 
decays £u£ with i = e, ^ are considered in both signal and background. 

The parton distribution functions used are MRS (A). 

Figure Captions 

[1] Feynman diagram topologies contributing at tree-level to the process bg — > 
M~i'~^i'iU£i, where i represents a lepton. Internal wavy lines represent a pho- 
ton, a Z or a W^, whereas the internal solid ones refer to a lepton, a neutrino, 
a bottom- or a top-quark, as appropriate. The total number of graph is 43 (ex- 
cluding Higgs couphngs). The single-top diagrams are number 2 and 3. Charge 
conjugated diagrams can be trivially obtained by reversing the fermion lines. 

[2] Differential distributions in (a) energy of the final state 6-jet (left plot) and its 
angular separation from the two leptons (right plot: from that generated by the 
top-decay W"^, solid line, and from that generated by the non-top-decay W"^, 
dashed line) and (b) transverse momentum (left plot) and pseudorapidity (right 
plot) of the final state 6-jet in events of the type bg — > tW^ — >• bW^W~ — >• 
b{l~V()[i'^V(ii), with £, £' = e, /X, at the LHC, (a) before the acceptance cuts and 
(b) after the acceptance cuts 1-6. The parton distribution functions used are 
MRS(A). 

[3] Differential distributions in the polar angle of the dilepton system with respect to 
the beam direction for the Higgs signal [Mh = 170 GeV, solid histogram) and 
the single-top background (dashed histogram) at the LHC after the acceptance 
cuts 1-6. The parton distribution functions used are MRS (A). Note that the 
background rates have been divided by two in order to facilitate the comparison 
between the two curves. 

[4] Differential distributions in the azimuthal angle of the dilepton system in the 
plane transverse to the beam direction for the Higgs signal {Mh = 170 GeV, solid 
histogram) and the single-top background (dashed histogram) at the LHC after 
the acceptance cuts 1-6. The parton distribution functions used are MRS (A). 
Note that the background rates have been divided by two in order to facilitate 
the comparison between the two curves. 

[5] Differential distributions in the estimated invariant mass of the i~ u^Uii sys- 
tem for the Higgs signal [M^ = 170 GeV, solid histogram) and the single-top 
background (dashed histogram) at the LHC after the acceptance cuts 1-6. The 
parton distribution functions used are MRS (A). Note that the background rates 
have been divided by two in order to facilitate the comparison between the two 
curves. 

[6] Differential distributions in the angle between the lepton with highest transverse 
momentum, boosted to the dilepton rest frame, and the momentum vector of the 



dilepton system itself for the Higgs signal {Mh — 170 GeV, solid histogram) and 

the single-top background (dashed histogram) at the LHC after the acceptance 
cuts 1-6. The parton distribution functions used are MRS (A). Note that the 
background rates have been divided by two in order to facilitate the comparison 
between the two curves. 
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hg tW^ bW+W- 
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13 


^ = 14 TeV 




£ = 5 fb-i 




mt = 


175 GeV 
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Diagrams by MadGraph 
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